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A simple model is studied for the chemotactic movement of biological cells in the presence of a 
periodic chemical wave. It incorporates the feature of adaptation that may play an important role 
in allowing for "rectified" chemotaxis: motion opposite the direction of wave propagation. The 
conditions under which such rectification occurs are elucidated in terms of the form and speed of 
the chemical wave, the velocity of chemotaxis, and the time scale for adaptation. An experimental 
test of the adaptation dynamics is proposed. 
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Many biological processes involve chemotaxis, cellu- 
lar motion in response to a chemical stimulus. Often, 
the chemo-attractant propagates through a set of cells 
as traveling waves [|l]j2| , as in a case of longstanding in- 
terest: the emergence of a multicellular structure from 
colonies of the eukaryotic microorganism Dictyostelium 
discoideum {Dd) In controlled experiments, a mono- 
layer with 10^ — 10^ cells/cm^ on the surface of agar 
begins within several hours after nutrient deprivation 
to support waves of cyclic adenosine monophosphate 
(cAMP) triggered by spontaneous release of cAMP from 
a small subpopulation of cells. These target or rotating 
spiral waves (Fig. 1), whose fronts appear as bands under 
dark-field visualization through their effects on cell shape 
Q, induce chemotaxis toward their centers, followed by 
complex multicellular morphogenesis. 

Chemical waves in excitable media such as Dd are quite 
thoroughly explored ||,^ , but their coupling to cell den- 
sity through chemotaxis is far less well-understood, al- 
though of longstanding interest iQ-^. As emphasized 
recently and illustrated in Fig. la chemotaxis 
driven by traveling waves is quite intriguing. A cell in the 
position indicated by the arrow experiences a progres- 
sion of leftward-moving wavefronts as the nearby spiral 
rotates outward. In seeking higher levels of cAMP, the 
cell would move first rightward into each advancing wave, 
then leftward after the peak has passed (Fig. lb). In the 
simplest model of chemotaxis, the cell velocity is pro- 
portional to the local chemical gradient, and it has been 
argued |10| (but not proven theoretically) that the net 
cellular motion would be in the same direction as the 
wave: i.e. "advection" away from the center, rather than 
the observed motion towards the spiral core. Tracking 
studies of cells ||l^ suggest a resolution to this by not- 
ing that as the cells experience the rising cAMP level of 
the approaching wave, their chemotactic response dimin- 
ishes, leaving them less responsive to the trailing edge, 
but their response recovers in time for the next front. 
They thus rectify the traveling waves, with net motion 
opposite that of the wave. 



In an effort to understand the underlying mechanism 
of this process, we study here a very simple model for 
"adaptive" traveling-wave chemotaxis and suggest ex- 
periments to test its predictions for the conditions un- 
der which rectification occurs. This model is closely re- 
lated to, but considerably simpler than one introduced 
recently in important work by Hofer, et al. |12|, who 
demonstrated by numerical computations that a process 
of adaptation could lead to rectified motion. A number 
of important aspects of this problem become clear with 
these simplifications, particularly in the experimentally 
relevant limit of chemotactic velocities small compared to 
the wave speed. First, in this limit an elementary proof 
is given of the heuristic argument ||l^ that nonadaptive 
chemotaxis will not produce rectified motion. Second, it 
is shown that rectified motion requires only two main 




FIG. 1. (a) Dark-field image of spiral waves in Dic- 
tyostelium discoideum Q]. Wavefronts of cAMP appear as 
dark bands. A cell in the position indicated experiences a pe- 
riodic train of cAMP waves, shown schematically in (b). Net 
chemotactic motion occurs toward the spiral core, opposite 
the direction of wave propagation. 
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ingredients: (i) a single characteristic time for adapta- 
tion, and (ii) a response function that decreases with 
concentration. 

Third, the net chemotactic flux is shown to have a ther- 
modynamic analogy in being proportional to the area en- 
closed by certain limit cycles exhibited in the response- 
concentration plane. Fourth, an analytical calculation 
confirms the intuitive notion that rectification is great- 
est when the adaptation time is comparable to the wave 
period (as seen in experiment pO[|), although there can 
be a delicate interplay between the competing processes 
of rectification and advection. Finally, an experimental 
test of these results is suggested. 

Consider a one-dimensional set of noninteracting cells 
at density p responding to a periodic chemical concen- 
tration wave c{x, t) = c{x -\- vt) with wavelength A and 
velocity v. Typically, A O.f — 0.5 cm (Fig. f), and 
w ~ f 0^ — f O'^ /im/min. We leave aside the complex dy- 
namics of wave production and its connection to the cell 
density [T^-[r^]. The model is formulated at the level of 
the coordinate x(t) of a cell, and for the present pur- 
poses is deterministic, as the random motions of the cells 
during one wave period are small on the scale of the wave- 
length A. Deterministic chemotaxis arising from chemical 
gradients is described by the overdamped dynamics 



dx d 

—— — r—c{x + vt) 
dt dx 



(1) 



The chemotactic response coefficient r measures the 
strength of chemotaxis, with cells migrating to high val- 
ues of c when r > 0. When r responds to c we have 
"adaptive chemotaxis" ; otherwise the motion is nonadap- 
tive. The flux of cells (J) (= (1/A) dzJ) averaged over 
one wave period is found by solving (|l]) in the moving 
frame z = x + vt, with dz /\r(dc/ dz) + v\ — dt, and trans- 
forming back (see also |l6|). 



J = {j)/pv = ([1 + r (dc/dz) /vr^y 
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(2) 



It is known from experiment that the typical chemotac- 
tic velocity r dc/dz in Dd is at least an order of magnitude 
lower than the wave speed v so we expand J in 

powers of v~^, 



^ I / dc\ 1 

dz 




(3) 



In nonadaptive chemotaxis, r is constant, so the first 
term vanishes by the periodicity of c{z). The first non- 
vanishing contribution to J7 is in the direction of the 
wave propagation, independent of the form of c(z) since 
the variance of rdc/dz is manifestly positive. This con- 
firms the heuristic argument of Wessels, et al. |l^]. The 
physical basis for this was emphasized in the context of 
Brownian particles forced by moving optical traps p6|: 



Particles migrating into the advancing wave experience 
the leading edge (and hence chemotax) for a shorter time 
than they do the trailing edge, since their velocity rela- 
tive to the wave is greater in the former case than in the 
latter. One step forward, two steps back, so to speak. 

A very simple adaptive chemotaxis model has two in- 
gredients: (i) an equilibrium "adaptation function" /(c) 
that is a decreasing function of c, and (ii) a single time 
constant r for the relaxation of r toward / [O , 



(4) 



A cell having experienced a concentration c for times 
much longer than r will have a response coefficient r — 
f{c) when next presented with a gradient: low when c is 
high, and vice versa. As c changes with time the response 
will attempt to equilibrate to /(c), but will lag behind 
when c changes on time scales shorter than t. One can 
think of this lag as a memory or inertial effect, and it 
provides a means of rectification. We expect the adap- 
tation time to be comparable to the refractory period of 
cAMP signaling. Eq. Q is one member of a FitzHugh- 
Nagumo model that has been studied in related work on 
chemotaxis p. 

Let us introduce the dimensionless coordinate Z = kz, 
time T — ujt, concentration C = cjc^, response coef- 
ficient i? — r/f^, and adaptation function F = f j f^., 
where fc = 2t:/\, lo = kv, cq is the peak wave concen- 
tration, and /o = /(c = 0). The rescaled dynamics are 



dz 

dT 
dR 
'dT 



SR 



dC 
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dZ 

= F{C) - R , 



with two dimensionless parameters, 

f2 — LOT , S — focok/v 



(5a) 
(5b) 

(6) 



The quantity Q measures the relaxation time in units of 
the wave period, while S is the ratio of a characteristic 
chemotactic speed to the wave speed. 

The expansion of the average flux in Eq. (H) now ap- 
pears as an expansion in S, with the 0{S) term possibly 
rectifying, and those of order always advective. For a 
given value of S, we expect three regimes of il: (i) fi ^ 1, 
the nonadaptive case already discussed, (ii) ~ 1, where 
rectification may occur, and (iii) 51 <C 1, with instanta- 
neous adaptation. In region (iii), the response tracks the 
concentration precisely, and the leading and trailing sides 
of the wave are not distinguished, so the fiux is positive 
and given analytically by Eq. (|^) with r replaced by 
/(c). It is smaller than in (i) since /(c) is smaller for 
high c. Rectification may occur in region (ii), where the 
down-regulation of the response triggered by the advanc- 
ing wave has not fully recovered by the time the trailing 
edge is encountered. 
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FIG. 2. Relation between morphogen concentration and 
chemotactic response. Panels show the imposed chemical 
wave C{Z{T)) (dashed) and the response coefficient R[Z(T)) 
(solid) for S = G.2 and various scaled frequencies fl. 

Figure ^ shows for different values of VL the concentra- 
tion C and response coefficient R as functions of time 
along the trajectory of a moving cell obtained by numer- 
ical integration of (|^). An arbitrary initial condition de- 
cays in a time of order 57 into the steady patterns shown. 
To simulate the sharply-peaked traveling waves seen in 
experiment, we chose C{Z) = exp [(5 (sin(Z) — 1.0)], with 
[3 = 3.0. The response function is the simplest: a linearly 
decreasing function of C: F{C) = 1 — C. For = 0.1 
(close to instantaneous adaptation) we see R almost pre- 
cisely anticorrelated with C, whereas for 17 ~ 1 the asym- 
metric response between leading and trailing edges is 
quite apparent. For = 10 the response settles to a 
nearly constant value determined by the mean value of 
the concentration over one wave period. 

The extent to which the response is "out of equilib- 
rium" is seen with limit cycles in the R — C phase plane 
shown in Fig. H, with position Z as a parameter. For 
small the cycle hugs the equilibrium curve F{C), while 
when f2 3> 1 it is a narrow loop encircling a horizontal 
line of constant R. In the rectifying regime (fi = 1) the 
cycle lies very far from equilibrium, forming a large closed 
loop. Using the high- velocity result in Eq. (||), we may 
express the chemotactic flux directly in terms of the area 
enclosed by this loop C, 



— / RdC , 



(7) 



in much the same way as we associate mechanical work 
with loops in the pressure-volume plane. With the sense 
of traversal of the loops shown in Fig. || this area is posi- 
tive, and hence rectifying. With these results, we obtain 
the chemotactic flux shown in Fig. ^, highlighting the 
window around 17 = 1 within which rectification occurs. 
For S small rectification occurs over almost the entire 
range of 17, since advection is negligible. But for larger 
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FIG. 3. Phase portrait of traveling-wave chemotaxis. Data 
in Fig. ^ are replotted in the R — C plane. The equilibrium 
function F{C) is shown dashed. Arrows indicate sense of 
traversal with increasing time. 

S advection dominates at the extremes of 17, and rectifi- 
cation occurs in a very narrow window near li = 1 . The 
qualitative behavior seen in Figs. is unchanged by 
the inclusion of more complicated nonlinear adaptation 
functions more faithful to the biochemistry of receptor 
binding fTl. 
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FIG. 4. Flux as a function of relaxation time and chemo- 
tactic velocity. At low S, rectification occurs over a broad 
range of Q,, while at higher S advection is dominant except 
over a narrow range of f2 ~ 1. 

Insight into the quantitative behavior of the model may 
be obtained by an analytic calculation in the limit of large 



wave speeds (small S) |12|. We assume an expansion of 
the position Z and response coefficient R in powers of S: 
Z{T) = Z^°\T) + SZW{T) + • • ■ and R{T) = i?(o)(r) + 
SR'^^\T) + At order S"" in this moving coordinate 
system we obtain dZ^°')/dT = 1, so Z(")(T) = T. At 
order S we obtain the equation of motion of a particle 
due to a time-dependent force 

dZW 



dT 



i?(°)(T) 



dC{T) 
dT 



(8a) 



i?W(T) 



dr'e-(^-^')/"F(C(T')) 



(8b) 
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We continue with the simple model F{C) = 1 — 
C. If C{Z) has the Fourier expansion C{Z) ^ C + 
S^^i [<Xn Gxp{inZ) + C.C.], then after the transients the 
modes of R are shifted in phase from the wave, 



n=l 



.„e'("^-^")+c.c. 

v/1 + {nny 



(9) 



The phase shifts On satisfy tan0„ = nQ, and thus are 
nearly zero for instantaneous adaptation and tend to tt/2 
as f2 — > cx), more rapidly with increasing mode number 
n. To compute the flux in this small-S* limit, we find the 
mean value of i?(o)(r)dC/dZ over one period, 



J 



S 

— y'nsin(26l„) \a„\ 



(10) 



where sin(20„) = 2nil/ (1 + (nfl)'^). The generic behavior 
of this sum is seen from its first term, which vanishes 
as J7 — > and — > oo, and peaks at 17 = 1 to give 
maximum forward flux. When the wave has dominant 
spectral weight in a higher mode n, the peak in flux will 
occur for ~ 1/n, as in Fig. 4. Once we go beyond the 
limit of vanishing S, the negative contributions to the 
flux (^ compete with the rectifying part to produce the 
resonance-like behavior seen in Fig. ^. 

We see that rectified chemotactic motion requires only 
two simple ingredients: a single time scale for adaptation, 
and an adaptation function that decreases with concen- 
tration. Adaptive phenomena are found in many bio- 
logical systems besides D. discoideum, including those 
of bacterial chemotaxis While this problem may 

appear similar to ones of directional transport studied 
recently in the context of electrophoresis , there is a 
fundamental distinction; the rectified motion arises here 
through processes internal to the particles, not through 
stochasticity. Indeed, the random motions of cells are ex- 
pected to decrease the efficiency of rectification. Finally, 
these results suggest experiments to probe the competi- 
tion between advective and adaptive chemotaxis (Fig. U) 
with artificially-produced chemical waves of controllable 
shape and velocity |^0| , complementary to recent experi- 
ments with fixed gradients pi]] and time- varying uniform 
concentrations |l^ . Observation of the predicted advec- 
tion with a wave whose frequency is either very small or 
very large compared to the natural signalling frequency 
would provide important evidence that adaptation asso- 
ciated with an internal time scale is the operative mech- 
anism in rectified chemotaxis. 

I am indebted to J.T. Bonner, E.G. Cox, K.J. Lee, R.H. 
Austin, and P. Nelson for many discussions, to P. Holmes, 
and C.H. Wiggins for suggestions, and to A. Goriely for 
pointing out the work of Hofer, et al.. This work was 
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